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Coronaviruses (CoVs) can cause trivial or fatal disease in humans and in animals. Detection methods for a 
wide range of CoVs are needed, to understand viral evolution, host range, transmission and maintenance 
in reservoirs. A new concept, “Multiprobe QPCR”, which uses a balanced mixture of competing discrete 
non- or moderately degenerated nuclease degradable (TaqMan®) probes was employed. It provides a 
broadly targeted and rational single tube real-time reverse transcription PCR (“NQPCR”) for the generic 
detection and discovery of CoV. Degenerate primers, previously published, and the new probes, were from 
a conserved stretch of open reading frame 1b, encoding the replicase. This multiprobe design reduced 
the degree of probe degeneration, which otherwise decreases the sensitivity, and allowed a preliminary 
classification of the amplified sequence directly from the QPCR trace. The split probe strategy allowed 
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TaqMan detection of down to 10 viral nucleic acid equivalents of CoV from all known CoV groups. Evaluation 
Avian coronavirus was with reference CoV strains, synthetic targets, human respiratory samples and avian fecal samples. 
Split probe strategy Infectious-Bronchitis-Virus (IBV)-related variants were found in 7 of 35 sample pools, from 100 wild 


mallards (Anas platyrhynchos). Ducks may spread and harbour CoVs. NQPCR can detect a wide range of 


CoVs, as illustrated for humans and ducks. 


© 2009 Elsevier B.V. All rights reserved. 


1. Introduction 


Coronaviruses (CoVs) are positive single-stranded RNA viruses. 
They belong to the family Coronaviridae in the order Nidovirales 
(Cavanagh, 1997). CoVs have the largest known nonsegmented viral 
RNA genome, 27-32 kb in length and is composed in its 5’-proximal 
two-thirds of two large open reading frames, ORFla and ORF1b, 
encoding the replicase complex (Brian and Baric, 2005). The CoV 
family traditionally is subdivided into three groups based on sero- 
logical and genetic properties (Cavanagh et al., 1993; Siddell et al., 
1983). Groups 1 and 2 infect a large range of mammalian species 
whereas group 3 is restricted to birds (Saif, 2004). Classification of 


Abbreviations: C,, threshold cycle; FRET, fluorescence resonance energy trans- 
fer; LNA, locked nucleic acid; NQPCR, nuclease-based real-time reverse transcription 
PCR; ORF, open reading frame; PCR, polymerase chain reaction; QPCR, quantita- 
tive real-time reverse transcription PCR; RNA, ribonucleic acid; RT-PCR, reverse 
transcriptase polymerase chain reaction; TET, tetrachloro-6-carboxyfluorescein; TD, 
touch down temperature protocol. 
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the severe acute respiratory syndrome coronavirus (SARS-CoV) in 
group 2 or as the prototype of a new group 4 is subject to con- 
troversy and is complicated by the seeming recombinant origin 
of its genome (Gorbalenya et al., 2004; Rest and Mindell, 2003; 
Snijder et al., 2003; Stavrinides and Guttman, 2004). The situation 
recently became more complicated because bats were found to har- 
bour a wide range of mammalian CoVs (Li et al., 2005; Poon et al., 
2005). It is thus likely that bats serve as reservoirs for CoVs in a 
similar manner as anseriform birds do for influenza A (Webster et 
al., 1992). CoVs are responsible for a broad spectrum of diseases, 
including respiratory and enteric pathologies, both in humans and 
in animals (Saif, 2004). Human CoVs(HCoVs) cause respiratory tract 
illnesses. However, there is also direct evidence for involvement 
in enteric (e.g. strain HECoV 4408) (Zhang et al., 1994; Zhu et al., 
2006) and suggestive evidence for involvement in neurological dis- 
ease (Clarke et al., 1979; Flewett et al., 1987; Foley and Leutenegger, 
2001; Gerna et al., 1985, 1984; Han et al., 2006; Jacomy et al., 2006; 
Luby et al., 1999; Resta et al., 1985; Saif, 2004; Schnagl et al., 1990, 
1986; Sitbon, 1985; Vabret et al., 2006; Wang et al., 2007; Zhu 
et al., 2006). HCoV are frequent causes of the common cold; usu- 
ally a self-limiting upper respiratory tract (URT) illness (Heikkinen 
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and Jarvinen, 2003). The prototypic HCoVs, HCoV-229E, and HCoV- 
OC43, were identified in the 1960s (Hamre and Procknow, 1966; 
Vabret et al., 2001). They are members of groups 1 and 2, respec- 
tively. Infection with these viruses is thought to be responsible for 
~30% of common cold cases (Fields et al., 2001). Until the discov- 
ery of the SARS-CoV, CoVs received relatively little attention as 
human pathogens beyond their causal role in the common cold. 
The most recently discovered HCoVs, the SARS-CoV (Peiris et al., 
2003), HCoV-NL63 (van der Hoek et al., 2006) and CoV HKU1 (Woo 
et al., 2005) causes severe respiratory tract infection. Addition- 
ally, the knowledge of the pathogenesis of CoVs in wildlife species 
is also limited. Among the many and highly diverse CoVs found 
in bats (see above), SARS-like CoVs were identified in horseshoe 
bats of the genus rhinolophus. Virological and serological studies 
indicated that masked palm civets (Paguma larvata) were infected 
with SARS-CoV or a closely related virus (Wang and Eaton, 2007). 
SARS-CoV thus may have spilled into civets from bats, and then 
emerged into humans. However, the lack of genomic sequences for 
many animal CoVs impeded clarification of the origin of SARS-CoV. 
Sequencing analysis of SARS-CoV revealed that the 5’ end, contain- 
ing the polymerase gene, is of mammalian origin, and the structural 
genes (excluding spike glycoprotein) at the 3’ end are of avian ori- 
gin (Jackwood, 2006). With some exceptions (Pantin-Jackwood et 
al., 2007; Spackman and Cameron, 1983; Spackman et al., 2005; 
Swayne et al., 2004; Woo et al., 2009), the spread of CoVs among 
birds is not well known (Cavanagh, 2005; Jonassen et al., 2005). 
Even though the number of known CoVs in birds is relatively small 
in comparison to the number in bats, like in bats, birds contain 
a wide diversity of CoVs. While this report was being prepared, 
three novel CoVs, classified as CoV group 3c, were identified in three 
different bird families, bulbul CoV HKU11, thrush CoV HKU11 and 
munia CoV HKU13 (BuCoV HKU11, ThCoV HKU12, MuCoV HKU13) 
(Woo et al., 2009). New diagnostic methods, which are able to detect 
a wide range of CoVs in various species and which can support the 
better understanding of CoV evolution, host range and reservoir 
conditions, are needed. 

Currently, diagnostic tests are mostly limited to sequence ampli- 
fication techniques as the CoVs have a restricted host range and are 
fastidious to propagate in cell culture. The clinical relevance of sero- 
logical methods, mainly based on enzyme-linked immunosorbent 
assays, is limited, since they are generally time-consuming. A vari- 
ety of assays for detection of CoV have been developed to obtain 
more sensitive and rapid diagnostic results, such as reverse tran- 
scriptase PCR (RT-PCR) (Vabret et al., 2001), nested RT-PCR (Myint 
et al., 1994) and, recently, also real-time RT-PCR (van Elden et al., 
2004). For simplicity, real-time RT-PCR is in the following referred 
to as “QPCR”. QPCR is less time-consuming than nested or unnested 
RT-PCR methods, because the amplification and analysis are com- 
pleted in a closed system, reducing contamination risks. 

Although molecular techniques enable the detection of species 
or groups of CoVs, there are few reports of pan-CoV RT-PCRs 


(Ksiazek et al., 2003; Moes et al., 2005; Stephensen et al., 1999), 
especially those that have been adapted to a real-time format 
(Escutenaire et al., 2007). All methods are based on amplification 
of fragments from the replicase gene, which has a highly con- 
served structure and function. This region is suitable for design of 
primers allowing a broad-spectrum detection of genetically distant 
CoVs. The SYBRgreen method (Escutenaire et al., 2007) is both sen- 
sitive and specific according to our experience. The specificity is 
ascertained by a melting curve analysis, which in many cases can 
discriminate a true positive from a false signal (Escutenaire et al., 
2007). However, a weakly positive true signal may occasionally be 
obscured by a stronger false one. Thus, SYBRgreen techniques to 
some degree lack the increased specificity provided by TaqMan® 
nuclease-based systems (Bustin, 2004). We therefore set out to 
develop a broadly targeted TaqMan® detection method system for 
Coronaviridae using a set of degenerate primers previously devel- 
oped and reported by us (Escutenaire et al., 2007) for the generic 
detection of CoVs. The availability of an alternative hybridization- 
based assay is useful for verification of results based on the previous 
method. The greater specificity may also facilitate the interpreta- 
tion of weak hybridization signals, which might be overlooked in 
the SYBRgreen technique. An optimised mixture of three probes, 
giving a signal arising from the best binding probes, was neces- 
sary to cover most of the known variation within the target stretch 
of Coronaviridae. They target discrete, short, conserved domains in 
coronaviral subsets. Their use avoided the extensive degeneration 
necessary with only one probe. LNA (locked nucleic acid) was used 
to increase the melting temperature of the relatively short probes. 
This enhances their affinity to the consensus target sequences. A 
somewhat similar approach has been used by others (Yip et al., 
2005), but without the broad, virus family-specific, targeting and 
the increased specificity provided by the nuclease-based probe 
detection of the present method. 

The goal of this study was to design a broadly targeted tech- 
nique for detection of all known groups of CoVs in animals and 
humans for use in a single tube QPCR. It is intended to complement 
or replace the previous SYBRgreen pan-CoV QPCR (Escutenaire et 
al., 2007) published by our group. The new, broadly targeted, PCR is 
in the following referred to as the nuclease-based QPCR (NQPCR). 
It provides a probe-based positive identification that increases the 
specificity of a positive result and a provisional classification when 
the three probes are tagged with different fluorophores. 


2. Materials and methods 


2.1. Computer program for identification of conserved portions of 
variable viral genomes 


The program, “ConSort”® (J. Blomberg et al., unpublished), 
uses alignments of 2-10000 viral sequences made using BLAST, 
CLUSTALW, MegAlign or MultAlin. Areas of sequence conserva- 


Fig. 1. (a) Alignment of the triple-probes with sequences from reference strains of CoV. Conserved regions are black. The most frequent variants are shown in gray. Variable 
positions have no background. The three probe sequences with IUPAC ambiguity codes are shown at the bottom, LNA positions are underlined. (b) Predicted hybridization 
of amplimer stretches from a broad range of coronaviruses to the three probes. An amplimer sequence from the Ottenby duck samples, with an IBV-like sequence (Ottenby 
67433, obtained from a mallard. It was not mentioned in Section 2. It was however identical to the Ottenby spring 3 sample in 120 of 130 positions) was used as query in a 
BLAST search for similar nucleotide sequences in GenBank. The figure lists the BLAST hits in order of decreasing similarity to the query. The stretches corresponding to the 
query found by the BLAST search were subjected to analysis with the Visual OMP (VOMP) software, which can predict the degree of interaction (shown as an absolute AG 
value) of an oligonucleotide with another. VOMP automatically removes all identical sequences, which meant that some reference sequences were removed from the hit list. 
For example, the OC43 sequence was removed. It was identical to the human enteric CoV sequence. Remaining reference sequences are shown in bold. All variants of each of 
the three probes were tested against all BLAST hit sequences. For each probe, the average absolute AG of the three variants predicted to hybridize best to the respective target 
sequences is shown. The three lines depict the predicted binding of probe_I (229E-derived; dark blue), probe_II (SARS-derived; red) and probe_III (NL63-derived, yellow), 
respectively. The columns show the average of the predicted probe absolute AG values, to indicate whether the triple-probe system would be able to detect that target 
sequence. (c) Differential fluorescence signal obtained with the three separately labeled probes. CoV, either RNA from reference strains, or synthetic target oligonucleotides 
from published CoV strains, were run in NQPCR. Results were colour coded as in (b). For further explanation of strain names, see Table 1. TCoV: Turkey CoV; CoV-cons-synth: 
synthetic oligonucleotide made from a consensus sequence of the PCR amplimer stretch (see Section 2). (For interpretation of the references to color in this figure legend, the 


reader is referred to the web version of the article.) 
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tion were identified, and the degree of representativeness for all 
sequences in the alignment computed. The program shows first the 
most common variant for each position, then less common ones in 
decreasing order. Insertions and deletions are marked especially. A 
graphical overview of the degree of conservation is also given. 


2.2. Primer and probe design and synthesis 


A set of previously reported primers were used (Escutenaire 
et al., 2007). Oligonucleotides were designed with the aim 
of (i) primers having an approximately equal Ty, (ii) the 
probe having a higher T,, than the primers, (iii) minimiz- 
ing primer-primer and primer-probe interactions, while (iv) 
still operating within a conserved region. All oligonucleotides 
were checked for homo- and heterospecific annealing as well 
as self-annealing loops, using a cutoff of —5kcal/mol for the 
predicted Gibbs free energy [AG] of oligonucleotide interac- 
tions. The primer-primer interactions were studied using Oligo 
Analyzer 3.0 (http://207.32.43.70/biotools/oligocalc/oligocalc.asp), 
and http://www.idtdna.com/analyzer/Applications/OligoAnalyzer. 
Primers were purchased from Thermo Hybaid, Interactiva 
Division (Ulm, Germany). The probes with LNA modified 
nucleotides were purchased from Eurogentec (Seraing, Belgium). 
The primer sequences and their positions in the genome of SARS- 
TOR2 (AY274119) are: 11-FW: 5’/-TGATGATGSNGTTGINTGYTAYAA- 
3’ (nt 15467-15670) and 13-RV: 5’-GCATWGTRTGYTGNGARCARAA- 
TTC-3’ (nt 15801-15825). The probe sequences and their position 
in the genome of HCoV-229E (GenBank accession No. NC0O02645, 
nt 14836-14862), SARS-TOR2 (GenBank accession No. AY274119, 
nt 15743-15770) and NL63 (GenBank accession No. NC005831, 
nt 14788-14804) are probe_I: 5’-[FAM|TTGTATTATCAGAATGGYGT- 
STTYATG[EDQ]-3’, probe_II: 5’-[ FAM |TGTGTTCATGTCWGARGCWA- 
AATGTTGG[EDQ]-3’, probe_III: 5’-[FAM|]TCTAARTGTTGGGTDGA- 
[EDQ]-3’. In later experiments, the same probe sequences were 
labeled with separate fluorophores. Then, probe_I was labeled with 
ROX and probe_III with JOE instead of FAM. The black hole quencher 


Table 1 
Collection of viral strains used for the evaluation of the pan-CoV NQPCR. 
Coronavirus species Strain Source 
Group 1 
HCoV-NL63 NL63 L.van der Hoek 
HCoV-229E 229E ATCC 
CCoV type | Elmo/02 C.Buonavoglia 
CCoV type Il 144/01 C.Buonavoglia 
FCoV type I FIPV DF2 I. Kiss 
PRCoV PRCV I.Kiss 
TGEV Purdue I.Kiss 
PEDV PEDV 7 I.Kiss 
Group 2 
HCoV-OC43 OC43 ATCC 
MHV MHV-A59 P.Rottier 
BCoV BCQ-4 A.Kheyar 
BCoV BCQ-3708 A.Kheyar 
Group 3 
IBV IBV 927 I.Kiss 
IBV 7/91 (793B) D.Cavanagh 
IBV HV-10 D.Cavanagh 
IBV Gray 390 D.Cavanagh 
IBV D207 D.Cavanagh 
IBV Arkansas 99 D.Cavanagh 
IBV UK/918/67 D.Cavanagh 
IBV AZRI 5508/95 D.Cavanagh 
IBV B1648 D.Cavanagh 
Group 4 
SARS-CoV FFM-ic M.Niedrig 


Reference 


EDQ was retained on these modified probes. Underlined positions 
denote LNA nucleotides. Introducing LNA residues at the most con- 
served positions was intended to increase the expected T;, of the 
probes, to achieve approximately the same Tj; for all three probes. 

The aim was to make the probes more resistant to viral varia- 
tion, allowing hybridization in spite of a few mismatches. However, 
probe-based techniques inherently have a greater specificity than 
non-probe-based (SYBRgreen) QPCRs. Thus, there is normally a 
trade-off between specificity and broadness of detection. In epi- 
demiological studies, with samples containing many different kinds 
of CoVs, several pan-CoV detection systems should be used to 
ensure a more complete coverage of the virus family. Alignment 
of the triple-probes with sequences from reference strains of CoV 
was generated using CLUSTALW version 1.83 (Fig. 1a). 

The predicted Gibbs free energy [AG] of oligonucleotide and 
target interactions was analysed utilizing Visual OMP'™ v6.6.0.0 
(DNA Software; http://www.DNAsoftware.com). The software uses 
nearest-neighbour thermodynamic theory to simulate hybridiza- 
tion of nucleic acids in solution and provides structural and 
thermodynamic parameters in intra- or intermolecular situations 
(SantaLucia, 2007). Data entry included the sequences of the 
three probes and target (the sequence stretch of CoV strains that 
hybridizes with the probe). The predicted free energy for hybridiza- 
tion of probes with different CoVs as target was obtained for all 
variants of each probe (i.e. the probes have different degrees of 
degeneracy; probe_I: 8, probe_II: 8 and probe_IlI: 6). The following 
solution conditions were used in Visual OMP to simulate the exper- 
iments; assay Tm 50°C, Na* 0.05, Mg2* 1.5 mM, 100 nM of probe and 
100 nM target. Under the above experimental conditions the aver- 
age AG for each probe variant with its target was plotted (Fig. 1b) 
(cf. its experimental counterpart, Fig. 1c). 


2.3. Collection of coronavirus strains 


A collection of CoVs was used to assess the sensitivity of the 
real-time RT-PCR (Table 1). Twenty-two human and animal strains, 


Copies/reaction Threshold cycle, G 


2.31 < 10° 12.96 
Hamre and Procknow (1966) 8.65 x 10? 22.14 
Pratelli et al. (2003a) 1.07 x 10° 14.23 
Pratelli et al. (2003b) 8.95 x 10® 10.74 
Evermann et al. (1981) 2.99 x 10° 12.54 

9.80 x 10! 29.50 
Bohl et al. (1972) 5.14 x 108 7.87 

7.94 x 104 18.50 
Kahn and McIntosh (2005) 3.03 x 10° 12,52 
Lavi et al. (1984) 2.46 x 107 9.08 
Milane et al. (1997) 9.30 x 10° 10.68 
Kourtesis et al. (2001) 4.21 x 10° 15.76 

4.42 x 10® ais 
Adzhar et al. (1995) 1.01 x 107 10.54 

1.14 x 107 1035 
Darbyshire et al. (1979) 6.01 x 10 26.52 
Davelaar et al. (1984) 7.51 x 10° 11.03 
Fields (1973) 1.33 x 108 13.87 
Cavanagh and Davis (1993) 4.67 x 10° 15.59 

1.73 x 107 9.65 
Shaw et al. (1996) 1.96 x 10° 17.02 
Drosten et al. (2003) 4.60 x 10! 30.75 


HCoV: human coronavirus; CCoV: canine coronavirus; FCoV: feline coronavirus; PRCoV: porcine respiratory coronavirus; TGEV: transmissible gastroenteritis virus; PEDV: 
porcine epidemic diarrhea virus; MHV: murine hepatitis virus; BCoV: bovine coronavirus; IBV: infectious bronchitis virus. 
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representative of the four CoV groups, including SARS-CoV, were 
analysed. The following references describe the origin of the collec- 
tion of CoV strains (Adzhar et al., 1995; Bohl et al., 1972; Cavanagh 
and Davis, 1993; Darbyshire et al., 1979; Davelaar et al., 1984; 
Drosten et al., 2003; Evermann et al., 1981; Fields, 1973; Hamre and 
Procknow, 1966; Kahn and McIntosh, 2005; Kourtesis et al., 2001; 
Lavi et al., 1984; Milane et al., 1997; Pratelli et al., 2003a,b; Shaw et 
al., 1996). 


2.4. Clinical samples and comparative methods 


Nasopharyngeal aspirates were collected from 77 patients 
admitted to the Uppsala Academic Hospital with an episode of 
acute respiratory infection. The samples were taken routinely for 
respiratory virus diagnosis by indirect immunofluorescence (IFA) 


Neighbor joining, bootstrap consensus 


of influenza A and B, parainfluenza 1, 2 and 3 and human res- 
piratory syncytial virus (the respiratory panel kit; cat. nr 3105, 
Chemicon International Inc., Temecula, CA, USA; or the Influenza 
A and B kits; cat. nr 5001 and 5002; Chemicon). For compari- 
son, a nested RT-PCR for influenza A and B was used as previously 
described (Herrmann et al., 2001). Nasopharyngeal secretion sam- 
ples were collected at the hospital by using a baby feeding tube 
and an aspiration trap. After suction the feeding tube was rinsed 
with approximately 2 mL of sterile saline and immediately trans- 
ported to the lab. Anjuman Hyder, Kare Bondeson and Bjorn 
Herrmann participated in the IFA diagnoses. Fecal samples of 100 
mallard ducks (Anas platyrhynchos) were collected (Ottenby Bird 
Observatory, Sweden) and pooled in 11 pools of autumn 2003 
and in 24 pools of spring 2004 and stored at —70°C until RNA 
extraction. 


Mallard Ottenby Autumn5 
Maillard Ottenby Autumn4 

Mallard Ottenby Autumn3 

Mallard Ottenby Autumn2 


Mallard Ottenby Spring13 


CoV group 3 


Mallard Ottenby Spring3 
Peafowl IBV AY641576 
Turkey CoV AF124991 
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62 | | Canine CoV AF 124986 


Human CoV HKU1 
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34]| Human Uppsala1903 
34" Bovine CoV mebus U00735 
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Fig. 2. Neighbour-joining bootstrap consensus tree, based on the short (179 bp) amplicon stretch. The bootstrap percentage of 500 trees is shown for most branches. The 
number of informative sites is low in this conserved portion of the RNA polymerase gene. Therefore, the branching pattern cannot be expected to precisely reflect the 
phylogeny and classification of Coronaviridae. The tree is intended to approximately represent the observed CoV sequences in a context of most similar reference sequences. 
The major CoV groups, and reference sequences belonging to them, are shown as boxes. Seven CoV positive human samples were sequenced (“Human Uppsala” plus sample 
number), and aligned with reference sequences. Four of them, which clustered with CoV group 2 strains, are shown in the tree. Sequences from three CoV positive samples 
were incomplete. Six of seven sequenced amplicons from mallards (“Mallard Ottenby spring/autumn” plus pool number) are also shown (see Suppl. Info, alignment in Fig. 
Sia and b, and two trees, Fig. S2a and b, made with other techniques, for further information). A more exact classification of the observed sequences will require sequencing 


of a longer stretch, which was out of scope for this methodological paper. 
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2.5. Sample processing and RNA extraction 


The samples were diluted 1:100 in nuclease-free water. RNA was 
extracted from 140 wl diluted sample using the QlAamp viral RNA 
kit (Qiagen, Hilden, Germany) according to the manufacturer’s rec- 
ommendations. RNA was recovered in 60 jl of nuclease-free water 
and either used immediately or stored at —70°C. 


2.6. Real-time QPCR system (“NQPCR”) 


A one-step real-time RT-PCR, nuclease-based QPCR (NQPCR), 
was developed based on the TaqMan® principle, using the iScript 
one-step RT-PCR kit for probes (Bio-Rad, Hercules, USA). The PCR 
was formed using a previously reported primer set (Escutenaire et 
al., 2007). All reactions were performed on the Corbett Research 
Rotor-Gene Real Time Amplification system (RG-3000, Corbett 
Research, Mortlake, NSW, Australia), with a total reaction vol- 
ume of 50 wl containing 25 pul 2x RT-PCR reaction mix for probes, 
800nM of each primer and 200nM of each probe, and 1 pl of 
iScript reverse transcriptase. The thermal conditions were as fol- 
lows; reverse transcription at 50°C for 30 min, followed by iScript 
reverse transcriptase inactivation and activation of the hot-start 
DNA polymerase at 95 °C for 10 min, followed by 5 touchdown steps, 
in which the annealing temperature was decreased from 56°C, by 
2 °Cevery second cycle down to 48 °C, followed by 50 cycles of 94°C 
for 30s and 46°C 60s. Fluorescence was measured during the latter 
period. The pan-CoV SYBRgreen QPCR (Escutenaire et al., 2007) was 
used as a control, both with the human respiratory, and the avian 
fecal, samples. 


2.7. Sequences of PCR products and phylogenetic analysis 


RT-PCR products from the SYBRgreen pan-CoV QPCR and NQPCR 
were gel purified (Qiaquick PCR purification kit; QIAGEN Inc.) 
and sequenced using the fluorescent dye terminator method, 
ABI PRISM® Big Dye!™ Terminator Cycle Sequencing v3.1 Ready 
Reaction kit (PerkinElmer) and on an ABI PRISM® 310 genetic ana- 
lyzer according to the manufacturer’s recommendations (Applied 
Biosystems, Foster City, CA, USA). The 179 nucleotide long amplicon 
sequences from mallard ducks from Ottenby and human nasopha- 
ryngeal aspirates from Uppsala were aligned using CLUSTALW 
Version 1.83 (Suppl. Info, Fig. Sla and b). Neighbour joining 
bootstrap consensus (Fig. 2) and minimum evolution, bootstrap 
consensus and maximum parsimony, bootstrap consensus (Suppl. 
Info, Fig. S2a and b) trees were created from the alignments using 
Mega 4.1 (Tamura et al., 2007). The reference coronaviral sequences 
were obtained from GenBank, with the given accession numbers. 


2.8. Quantitation standards 


Standard curves were generated for quantitation of assays with 
the help of three kinds of standard: (i) bovine CoV RNA, (ii) bovine 
CoV cDNA, and (ili) a synthetic CoV consensus oligonucleotide. 
The bovine CoV RNA is described in Table 1. The bovine cDNA 
was obtained by reverse transcription from this source. The RNA 
and DNA concentrations were determined spectrophotometrically 
using the Nanodrop ND_1000 instrument according to the manu- 
facturer’s instructions (Nanodrop, Wilmington, DE). 

Standard curves were also generated using a 10-fold dilution 
series of a synthetic CoV consensus target DNA oligonucleotide 
(“CORONA_CONSENSUS”: TGATGATGGTGTTGTGTGTTATAACAATGT- 
TIGTATTATCAGAATAATGTGTTCATGTCIGAAGCIAAATGTIGGGTTGA- 
ACCAGACATAAATAAAGGACCTCATGAATTTTGTTCACAGCATACAAT- 
GC), where the first underlined sequence corresponds to probe_I, 
the italicised to probe_II and the second underlined sequence to 
probe_III. Standards for standard curves were diluted in RNAse 


free water including 20 ng/wl of yeast tRNA (Ambion, Huntingdon, 
United Kingdom). 

The concentration interpolated from the standard curves was 
referred to as “equivalents” per PCR reaction. 


3. Results 


3.1. Identification of conserved portions in coronaviral genomes 
by using the ConSort® program 


A ConSort® analysis of a CLUSTALW alignment of 32 full-length 
coronaviral sequences indicated that ORF1b was one of the most 
conserved portions. During the design of the NQPCR probes, repre- 
sentatives from three CoV groups were used as query in BLAST at 
the NCBI website to search and align the viral sequences in Gen- 
Bank, i.e. probe_I: ORF1b of 229E; probe_II: ORF1b of SARS-TOR2; 
probe_III: ORF1b of NL63 were used as query. The ConSort output 
from the GenBank search was used for design of probes for QPCR. 


3.2. Evaluation of the NQPCR on coronavirus strains 


A major problem with nuclease-based TaqMan® QPCR probes 
is the sensitivity to mismatches. A in non-degenerated probe (5- 
[FAM ]TATTATCAAAATAATGTCTTTATG|TAMRAJ-3) was first synthe- 
sized. Not unexpectedly, it did not detect all CoV strains. A degen- 
erated probe 5’-[FAM|TAYTAYCARAAYRRYGTBTTYATG[|[TAMRA|]-3’, 
encompassing the whole variation, was then synthesized, with a 
degeneracy of 768. In our experience, probes with a high degener- 
acy give very low signals. In order to bring down the degeneracy per 
probe, three probes were designed. The probes were derived from 
three CoV groups, as described in Section 2. Probe_I and probe_lll 
were overlapping. The probes had a degeneracy of 8, 8 and 6, respec- 
tively, and LNA nucleotides were introduced in order to increase the 
Tm of two of the probes, to achieve approximately same Tj, for all 
three probes. When they were tested singly, they gave high relative 
signal strength, around 0.5, and detected most, but not all (TGEV, 
FCoV type I, IBV, CCoV type I, PRCoV, SARS, CoV-NL63, HCoV-OC43, 
BcoV11 and BcoV22), coronaviral strains. However, a NQPCR with 
a 1+1+1 mixture with 200nM of each FAM-labeled probes_l, II 
and III gave somewhat lower signal strength (0.3), but detected all 
coronaviral strains. However, when the three probes were labeled 
with separate fluorophores, signal strengths varied from 0.25 up to 
1. This probe mixture was predicted to bind to most, possibly all, 
known CoV strains (Fig. 1a and b) and was used in the further exper- 
iments. A varying pattern of predicted binding was obtained when 
reference and divergent target CoV sequences were tested in VOMP. 
No strain came out completely negative in the analysis (Fig. 1b). 
The predicted hybridization values were largely corroborated when 
strains were tested experimentally (Fig. 1c). 

Twenty-two human and animal strains representative of the 
four CoV groups, including SARS-CoV, were analysed and detected. 
The spectrum of detection of the assay was further verified by 
amplifying cDNA synthesized from RNA extracted from 10 CoV 
strains; CCoV, PRCoV, HCoV-229E, IBV 297, FCoV, TGEV, Bcov-22, 
SARS, IBV, BCoV-11, and HCoV-NL63. All of them were detected by 
the three-probe NQPCR system (Table 1). 


3.3. Analysis of human respiratory samples 


RNA from 77 human nasopharyngeal aspirates was analysed 
by both the NQPCR and the SYBRgreen QPCRs. Both techniques 
indicated that eight samples contained CoV. The SYBRgreen QPCR 
products of the eight CoV positive samples were sequenced, and 
aligned with reference sequences. Four of the seven sequenced 
amplicons turned out to belong to CoV group 2, being OC43/MHV- 
like (Fig. 2). One belonged to CoV group 1 and was HCoV-NL63-like 
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(Suppl. Info, Fig. S1a). However, it, and two more was incompletely 
sequenced. They were not included in the trees. 

No signal was obtained in the following RNA controls; 3 
enterovirus, 2 calicivirus and 2 influenza A positive samples estab- 
lished in previous work. Another sign of specificity was that among 
the 77 respiratory samples there were 5 RSV positive (by IFA) sam- 
ples. All 5 were negative in NQPCR. 


3.4. Analysis of avian fecal samples 


RNA from 35 pooled duck fecal samples in total representing 
100 wild birds, 24 pools from spring and 11 from autumn was anal- 
ysed by both the NQPCR and the SYBRgreen QPCRs. Of the 35 pools, 
5/11 autumn pools (Mallard_Ottenby_Autum 2, 3, 4,5 and 11) and 
2/11 spring pools (Mallard_Ottenby_Spring 3 and 13) were CoV pos- 
itive by NQPCR. Two of the CoV positive pools (one of the spring 
and one of the autumn pools) were also positive by SYBRgreen 
QPCR, but five of the seven CoV positive pools were missed by the 
SYBRgreen QPCR. The PCR products of the seven CoV positive sam- 
ples were sequenced (from both forward and reverse strands). The 
sequencing of one pool came out incomplete. All six fully sequenced 
amplicons were aligned with reference sequences. They resembled 
the peafowl CoV, accession number AY641576 and avian infec- 
tious bronchitis virus, accession number NC_001451, with sequence 
identities ranging from 83 to 85% and 84 to 86% (Suppl. Info, Fig. 
S1b). The sequences clustered together, indicating the existence of 
new CoVs in ducks. The seventh, incompletely sequenced, sample 
yielded a similar sequence. A caveat is the conserved nature of the 
targeted stretch. It is not suitable for fine discrimination of coronavi- 
ral strains. The major groups do however cluster approximately as 
expected. Neighbour joining bootstrap consensus (Fig. 2) and min- 
imum evolution, bootstrap consensus and maximum parsimony, 
bootstrap consensus (Suppl. Info, Fig. S2a and b) trees were created 
from the alignment. The branching pattern was almost identical in 
all three trees. However, further sequencing of longer stretches is 
necessary to establish the duck corona sequences as “new” avian 
coronaviruses, cf. (Woo et al., 2009). Sequencing of a longer stretch 
(559 nt without primer sequences, which include the stretch used 
by Woo et al. (2009) of the polymerase gene of one of the isolates 
was performed (Mallard_Ottenby_Spring 13). A comparison of the 
percent identity of Mallard_Ottenby_Spring 13 nucleotide sequence 
with published sequences of IBV strains and other CoVs is presented 
in Suppl. Info, Fig. S3. The similarity matrix revealed a percentage 
identity of Mallard_Ottenby_Spring 13 to IBV strains ranging from 
72 to 74%. The IBV strains possess a nucleotide identity of 88 to 92% 
within the group 3 (IBV proper) of CoVs, and nucleotide identity 
60% to other CoVs. This implies that the Mallard_Ottenby_Spring 13 
is an IBV-like sequence which is relatively dissimilar from estab- 
lished IBV strains. A more comprehensive analysis will be given in 
a forthcoming manuscript. 


3.5. Calculation of specificity and sensitivity of NQPCR 


The molecular limit of detection of NQPCR was 1-10 DNA copies 
of synthetic DNA nucleotides from a consensus CoV sequence. 
Serial 10-fold dilutions, from 10° to 10°, of the synthetic target 
was prepared in RNAse free water including 20 ng/wl of yeast tRNA 
(Ambion, Huntingdon, United Kingdom) as carrier, to determine the 
lower limit of detection of the assay. The amplification efficiency of 
a 10-fold serial dilution of BCoV RNA sample with an unknown con- 
centration down to negativity was also determined. The standard 
curves displayed a linear relationship for NQPCR for the assay, cor- 
relation coefficient (R2) of 0.99 and efficiency of 93%, and an R2 of 
0.99 and efficiency of 83% (Fig. 3a and b). Thus, the CoV DNA and 
CoV RNA targets gave similar efficiencies. 
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Fig. 3. (a) Amplification plot and standard curve, describing sensitivity tests. Primers 
11-FW/13-RV were used in combination with FAM-labeled probes_l, II and III. The 
curves represent, from left to right, 10’-10° copies of synthetic corona consensus 
130 oligonucleotide, with an efficiency of 0.93 and R? 0.99. No signals were obtained 
in the negative control. (b) Amplification plot and standard curve, describing an 
amplification efficiency test. Primers 11-FW/13-RV were used in combination with 
FAM-labeled probes_l, IJ and III, and run on 10-fold serial dilutions of BCov RNA, with 
an efficiency of 0.83 and R? =0.99. Thus, amplifying from CoV RNA gave a similar, 
but somewhat lower, efficiency. 


The sensitivity of the NQPCR was studied further by end-point 
titration of cDNA from BCoV with 6 observations of each dilution 
step, into the stochastic zone of 1-10 copies per NQPCR reaction. 
The cDNA sample was titrated in 10-fold dilution steps down to 
complete negativity. This sensitivity test result indicated a sensi- 
tivity of 1-10 copies per PCR assay, with a clear stochastic zone of 
positivity (Table 2). 

The sensitivity of NQPCR was the same regardless of whether 
separately labeled or commonly FAM-labeled probes were used 
(data not shown). 


Table 2 
NQPCR quantification result of 10-fold serially diluted BCoV RNA into the stochastic 
zone (1-10 copies per PCR reaction). 


Estimated copy number per PCR reaction 


100 10 1? 0 
Dilution 1 107 0 0 0 
Dilution 2 95 10 0 0 
Dilution 3 119 0 0 0 
Dilution 4 67 16 0 0 
Dilution 5 63 0 0 0 
Dilution 6 109 0 13 0 
S.D.> 23.27 6.98 5.31 0 


Six different tubes were analysed per dilution. The estimated RNA concentration in 
copies per PCR reaction, in relation to the Threshold cycle (C,) is shown. 

¢ At the estimated concentration of one copy per PCR reaction one of six samples 
were positive. 

> Standard deviation. 
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3.6. Reproducibility of the assay 


The reproducibility of the NQPCR was evaluated both inter-assay 
and intra-assay. The intra-assay variability was assessed using high 
and low amounts of cDNA from 10 different CoV strains. The 10 sam- 
ples were run in three replicates in same NQPCR run; the intra-assay 
variability of number of equivalents derived from interpolation in 
the standard curve had a standard deviation of 0.54 and a standard 
error of the mean (S.E.M.) of 10%. Inter-assay variability was deter- 
mined by testing the same set of cDNA samples in three different 
runs, giving a standard deviation of 0.45 and a S.E.M. of 8%. Judging 
from these results, at least 1 C; (two standard deviations) is needed 
to differentiate between the values of two samples. 


3.7. RNA extraction and RNA integrity controls 


As control for errors occurring after the receipt of samples (post- 
receipt errors), 11 CoV negative samples from ducks (belonging to 
the same series as the CoV positive samples) were spiked with a 
known concentration of IBV. They gave a CG value of 8.0 S.D. 0.3. 
When the same amount of IBV sample was added to three samples 
of phosphate buffered saline, a signal with the same G value 8.0 
S.D. 0.1 was obtained. Thus, there was no detectable inhibition of 
the PCR in the CoV negative samples. 

To monitor pre-receipt sample quality, all 35 duck fecal pool 
samples from Ottenby samples were analysed utilizing a histone 3.3 
PCR, which is based on evolutionarily conserved target sequences 
in histone 3.3, and can be run on samples from most vertebrates, 
including birds (Andersson et al., 2005; Forsman et al., 2005; 
Muradrasoli et al., 2006). The histone 3.3 PCR was run with (+rt) 
and without (—rt) the presence of reverse transcriptase. The —rt 
results signify the presence of eukaryotic (presumably from duck) 
DNA, whereas the difference between +rt and —rt results signifies 
the presence of eukaryotic RNA. The —rt reaction frequently came 
out negative, while the +rt reaction was positive in most samples. 
This means that histone 3.3 RNA was more abundant than histone 
3.3 DNA, which is logical. RNA often occurs in great excess over DNA 
in eukaryotic cells. The means and standard deviation (mean + S.D.) 
of the C,-value of histone 3.3 PCRs were calculated. The 28 CoV 
negative, histone 3.3 RNA positive samples yielded an average C,- 
value of 37.5 + 1.6, whereas the seven CoV positive and histone 3.3 
RNA positive samples yielded an average C;-value of 37.2 + 0.5. Stu- 
dent’s t-test did not detect significantly deviating histone 3.3 +rt 
and —rt CG values between CoV positive and CoV negative samples 
(p-value = 0.758). Moreover, one of seven CoV positive samples, and 
four of 28 CoV negative samples were histone 3.3 +rt negative. The 
difference between 1/7 vs. 4/28 is not significant (Fisher’s exact 
test; 2-tailed; p-value=1). Thus, there were approximately equal 
amounts of eukaryotic RNA in the CoV positive and CoV negative 
duck fecal samples, making it unlikely that the CoV negativity of 
some samples was due to RNA degradation. 


3.8. Comparison of touchdown and non-touchdown NQPCR 


Formation of non-specific products might interfere with ampli- 
fication of the specific target in low-copy-number samples. An 
advantage of the presented NQPCR is that the initial high anneal- 
ing temperature reduces amplification of non-specific PCR, and 
increases the specificity and sensitivity of the assay (Hecker and 
Roux, 1996; Mohamed et al., 2006). In developing the NQPCR assay, 
a standard non-touchdown PCR assay was first used. However, 
to further increase the sensitivity for aberrant targets a touch- 
down procedure was designed and assessed side by side with 
the non-touchdown PCR, by amplifying cDNA synthesized from 
RNA extracted from 10 CoV strains; CCoV, PRCoV, HCoV-229E, IBV 
297, FCoV, TGEV, BCoV-22, HCoV-229E, SARS, IBV, BCoV-11, and 


HCoV-NL63 (Bustin, 2004; Mohamed et al., 2006). The amplifi- 
cation efficiency was increased with the touch down procedure 
compared to the non-touchdown procedure. Non-touchdown had 
an efficiency of 70% (R2 0.99) while touchdown gave an effi- 
ciency of 90% (R2 0.99). Gel electrophoresis analysis of the PCR 
product indicated a higher yield compared to the non-touchdown 
procedure. 


4. Discussion 


CoV came into the searchlight when the SARS pandemic 
occurred in 2003. The intense researched spurred then has revealed 
an unexpected variety of CoV in various host species. However, 
many questions are still unanswered in CoV biology, including 
important aspects of viral evolution, subpopulations, hosts, vari- 
ation in pathogenicity and tissue tropism and regarding reservoirs, 
which are important in viral “survival” and transmission in nature. 
Bats are now recognised as reservoirs of a broad range of mam- 
malian CoVs (Poon et al., 2005; Wang and Eaton, 2007). Recently, 
several findings of CoV in wild birds have been reported (Barr et 
al., 1988; Jonassen et al., 2005; Liu et al., 2005). The results of this 
report extend these observations. They indicate that birds, another 
flying vertebrate, can act as reservoirs for avian CoVs. Considering 
the many unanswered questions regarding this and other aspects 
of coronaviral epidemiology, there is a need for novel methods, 
which are broadly detecting the various CoV variants in different 
hosts and which are also sensitive, simple and specific. Setting 
up a generic PCR requires an extensive analysis of variability, a 
time-consuming evaluation of candidate degenerate primers and 
probes, and a comprehensive collection of samples. In this study, a 
one-step real-time RT-PCR, was developed targeting a conserved 
region of the CoV ORF1b used previously in a SYBRgreen QPCR 
(Escutenaire et al., 2007). The Consort® program (Blomberg et al., 
to be published) was used. The new technique was based on nucle- 
ase probe degradation (TaqMan®) chemistry, degenerate primers 
and a low annealing temperature. The assay enabled detection of 
22 CoV reference strains, covering a large portion of Coronaviri- 
dae (Table 1). The coverage of the new bat CoVs is demonstrated 
by using synthetic target sequences of three most deviating HKU- 
CoV: HKU3-CoV, HKU5-CoV and HKU9-CoV) (Fig. 1c). They were all 
detected by the probe combination. 

Non-coronaviral controls were uniformly negative. A major 
problem with TaqMan® QPCR probes is the sensitivity to mis- 
matches. A nondegenerate probe based on the consensus of the 
intervening sequence between the primers was first synthesized. 
It did not detect all coronaviral strains. A degenerate probe encom- 
passing the whole variation, with a degeneracy of 768, was then 
tested. Not unexpectedly, it gave a very low signal. In our experi- 
ence (contrary to primers which can work even with a degeneracy 
of over 100) probes with degeneracy above 10 rarely give sufficient 
signal strength. In order to bring the number of variants down, the 
CoV target variation was analysed. It was found possible to spread 
it onto three overlapping TaqMan® QPCR probes of much lower 
degeneracy, having 8, 8 and 6 variants each. LNA residues were inte- 
grated in two of the probes (Mohrle et al., 2005; Nielsen et al., 2004), 
to enhance the Tj. The introduction of LNA allowed shortening of 
them, to focus on a short, relatively conserved, stretch (probe_l: 
27 nt, 8 LNA, probe_III: 17 nt, 6 LNA). Probe_II could be 28 nt long 
despite a low degeneracy. As demonstrated here, the probe mix- 
ture had a remarkable tolerance to mismatches. The free energy for 
hybridization was predicted using Visual OMP™. It demonstrated 
that hybridization between the three probes and a broad range of 
CoV targets was likely. This was further verified by experimental 
data that agreed with the result obtained using the Visual OMP™ 
(Fig. 1b and c). Interestingly, competition between the probes was 
not an appreciable problem. Although several of the probes some- 
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times gave a signal, the pattern of positivity indicated which group 
of CoV that gave rise to the signal. The mismatch tolerance was 
tested by synthetic target sequences with different degree of mis- 
match to the probes. Signals from synthetic DNA targets with 1, 2 
and 5 mismatches were easily detectable (data not shown). The 
probes were tested both with the same fluorophore (FAM), and 
with separate fluorophores. The sensitivity was the same for these 
alternatives. Although based on a short stretch in the conserved 
RNA polymerase gene, the latter alternative provided a preliminary 
identification of the CoV. 

PCR efficiency can be affected by numerous factors, inhibition 
in the sample, secondary structure interference, how well primers 
and probes fit to their targets, and how optimal the PCR condi- 
tions are. It is acknowledged that the efficiency calculation from 
the slope of a calibration curve often overestimates the PCR effi- 
ciency. We therefore determined the amplification efficiency of the 
NQPCR both for a 10-fold serial dilution of synthetic target and a 
10-fold serial dilution of BCoV RNA samples with an unknown con- 
centration down to negativity. Both standard curves displayed a 
linear relationship for NQPCR, with a correlation coefficient (R2) of 
0.99 and efficiency of 93%, and an R2 of 0.99 and efficiency of 83%, 
for the synthetic and natural, imperfectly fitting (BCoV), targets, 
respectively. This calls for caution regarding quantitative interpre- 
tation of C; values. Differences in efficiency of the PCR could also 
be due to the fact that RNA-to-cDNA conversion is dependent on 
template abundance. The possibility that mishandling of samples 
(RNA losses; pre-receipt errors) and PCR inhibition (post-receipt 
errors) contributed significantly to the results reported here was 
made unlikely by two control experiments. Signals from RNA of 
a highly conserved host housekeeping gene (histone 3.3) were the 
Same in CoV positive and CoV negative duck samples. This is a novel 
application of a control PCR which was earlier developed in the 
group (Andersson et al., 2005; Forsman et al., 2005; Muradrasoli 
et al., 2006). The difference between +rt and —rt signal levels from 
duck fecal samples were relatively even, indicating a rather constant 
admixture of eukaryotic RNA in them. 

The detection of HCoV-OC43-like and HCoV-NL63 sequences in 
8 out of 77 nasopharyngeal specimens of human origin, and duck 
CoV sequences among seven of the pooled wild bird fecal samples, 
provides evidence of the efficiency of the technique with clinical 
and field material. The pan-CoV SYBRgreen QPCR missed five duck 
samples that the NQPCR detected, since the cutoff (threshold) was 
set higher because of primer dimer signal. Primer dimers do not give 
a signal in a nuclease and probe-based (TaqMan®) technique like 
NQPCR. The NQPCR results from the human nasopharyngeal sam- 
ples were concordant with those of the pan-CoV SYBRgreen QPCR 
technique used as a control. The 8 NQPCR CoV positive samples 
were negative in the respiratory virus IFA, and in a nested influenza 
A+B PCR (data courtesy of dr BjOrn Herrmann), both of which do 
not detect CoVs. The method thus should be suitable for diagnostic 
purposes. 

Stephensen et al. (1999) described initially a set of consen- 
sus primers targeting a 251-bp fragment of the ORF1b. These 
primer sequences were subsequently modified for better reactivity, 
notably to the then newly identified HCoOV-NL63 (Moes et al., 2005). 
Another generic coronaviral PCR was recently published (Vijgen et 
al., 2008). Furthermore, (Sampath et al., 2005) reported a recently a 
broad-range PCR targeting the same conserved region of the ORF1b 
followed by electrospray ionisation mass spectrometry and base 
composition analysis for viral identification. Viral diagnostic tech- 
niques are undergoing a transition. Immunoassays are still valid 
for some purposes, but real-time PCR and microarray techniques 
are becoming more important. The technique presented above is 
simple, sensitive and is able to detect a wide variety of CoV. The 
extension of the assay by using three different fluorophores, with 
one for each of them, is a valuable modification. As expected from 


the alignment and predicted AG values, probe 3 was primarily 
active in positive avian (IBV-like) samples while probe 2 was posi- 
tive ina SARS sample (data not shown). Although not studied, probe 
3 should primarily detect CoV groups 1 and 2. 


5. Conclusions 


In conclusion, a QPCR with a novel strategy, which uses a mix- 
ture of LNA and non-LNA containing probes with an optimised low 
number of degenerations, is reported here. To the best of our knowl- 
edge, this is the first report describing such a strategy to reduce 
probe degeneracy while maintaining a broad sequence targeting. 
It should have a wide applicability. High sensitivity was indicated 
by titrations of a CoV consensus synthetic target DNA content and 
viral cDNA. The new NQPCR assay, positive with 22 highly diverse 
Strains of CoV, is a useful tool for demonstrating the role of human 
and wild bird CoVs in infections of the respiratory and digestive 
tract. This NQPCR thus proved its suitability also for the discovery 
of new, or so far unknown variants of CoVs. 
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